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ABSTRACT 

There is good evidence that most of the baryons in the Universe are dark 
and some evidence that most of the matter in the Universe is nonbaryonic 
with cold dark matter (cdm) being a promising possibility. We discuss ex- 
pectations for the abundance of baryons and cdm in the halo of our galaxy 
and locally. We show that in plausible cdm models the local density of cdm 
is at least 10~^^gcm~^. We also discuss what one can learn about the the 
local cdm density from microlensing of stars in the LMC by dark stars in 
the halo and, based upon a suite of reasonable two-component halo models, 
conclude that microlensing is not a sensitive probe of the local cdm density. 
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Introduction. While the quantity and composition of matter in the Uni- 
verse is still not known with certainty, it is known that: (i) luminous matter 
(stars, etc.) contributes much less than 1% of critical density; (ii) based upon 
primordial nucleosynthesis baryons contribute between about 1% and 10% 
of critical density; and (iii) based upon numerous dynamical measurements 
the total mass density is at least 10% of critical density [|l|, with several de- 
terminations indicating that it is close to the critical density |0 . Thus, there 
is overwhelming evidence that most of the matter in the Universe is "dark," 
compelling evidence that most of the baryons are dark, and mounting evi- 
dence that most of the matter is nonbaryonic [^. If the mean mass density 
is greater than about 10% of critical there are two dark matter problems, the 
nature of the baryonic and nonbaryonic dark matter. 

The case for a critical Universe with nonbaryonic dark matter receives 
further support from studies of structure formation: The most successful 
models rely upon nonbaryonic dark matter, and the cold-dark matter mod- 
els of structure formation (inflation-produced density perturbations and non- 
baryonic dark matter with negligible velocity dispersion) are very attractive 
[|]. The best motivated cold dark matter candidates are an axion of mass 
around 10~^ eV and a neutralino of mass between 10 GeV and 1 TeV |p, and 
large-scale experiments are underway to directly detect the axions or neu- 
tralinos in the halo 0. Needless to say, theoretical expectations for, and 
observational information about, the local mass density of cold dark matter 
are of great importance. 

The fiat rotation curves of spiral galaxies indicate that the luminous, disk 
shaped portion of a typical spiral sits in a dark halo that is roughly spherical 
with density that decreases as l/(r^ -|- a^) and extent that is undetermined 
(r is distance from the center of the galaxy and a is the core radius). The 
fact that galactic halos are more spherical and extended than the luminous 
parts of spirals strongly suggests that the dark halo material has probably 
not undergone significant dissipation. For the Milky Way, galactic modeling 
indicates that the core radius is between 2 kpc and 8 kpc and that the halo 
density nearby (r = ro — 8.5 kpc) is about 5 x 10~^^ g cm~^ (to within a factor 
of two) 0. At our position, the halo material supports around 130kms~^ 
of the 220kms~^ circular rotation velocity (the various contributions — halo, 
disk, etc. — to the rotation velocity add in quadrature.) 

In a cold dark matter Universe there are at least two forms of dark matter, 
baryons and cdm particles, and both are expected to contribute to the halo 
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mass density. The concerns of this Letter are twofold: first, the theoretical 
expectations for the local cdm density, and second, what one can learn about 
the local cdm density from microlensing experiments which can probe the 
baryonic component that exists in dark stars of mass 10~^Mq to IOOMq 
[§]. Since the EROS and MACHO microlensing searches now have candidate 
microlensing events 0, this is a very timely issue. 

Expectations. Consider a cold dark matter universe with ficdm + = 1- 
Based upon primordial nucleosynthesis O.Ol/i"^ < Qb < 0.02/i~^; further, 
in a critical, matter-dominated Universe the Hubble constant must be near 
its lower extreme, h ~ 0.5, in order to accommodate a sufficiently aged 
Universe (Hq = 100/ikms~^ Mpc~^ is the present Hubble constant and fij is 
the fraction of critical density contributed by species i). This means that the 
universal baryon fraction of the matter density is fs = = 0.04 — 0.1 and 
/cdm = 1 - ~ 0.9 - 0.96. 

(In the two popular variants of cold dark matter, hot + cold dark matter, 
where fl^ ~ 0.3, and cold dark matter + cosmological constant, where fl\ ~ 
0.8, flcdm ~ 0.2 and h ~ 0.7, the situation is only quantitatively different. 
In the former case, phase-space considerations limit the neutrino fraction 
in galaxies like our own to be less than about 5% and so the baryonic 
fraction of matter that can clump in galaxies, J'b — ^b/ i^—^iy) — 0.06—0.14, 
is slightly higher. In the latter case, the baryonic fraction of matter, Jb = 
0^/(1 — JIa) 0.1 — 0.2, is even larger. All that follows is applicable to these 
models by taking account of the larger value of /s.) 

So long as gravity alone shapes the evolution of the Universe the bary- 
onic fraction of matter remains at its universal value. Once dissipative forces 
(e.g., electromagnetic interactions) become important, baryons can lose en- 
ergy and become more condensed, increasing the local baryonic fraction. The 
formation of the disks of spiral galaxies provides a good example: Through 
collisional processes baryons lose their energy, but not angular momentum, 
ultimately forming a thin, rotationally supported disk. The local mass den- 
sity of the disk is about a factor of 20 greater than that of the halo. 

In turn we now consider three models for the formation of the halo of 
our galaxy, from a very simple scenario where only gravity is involved to an 
extreme scenario where hydrodynamical forces form the galaxy and the cold 
dark matter particles are captured subsequently by accretion. 

(1) The simplest and most plausible scenario is one where the halo of 
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our galaxy formed through the action of gravity alone, after which a small 
fraction of the baryons dissipated energy forming the disk. One then expects 
an isothermal halo, with density decreasing as r^^, whose extent is dictated 
by the "sphere of influence" of our galaxy, at most about half the distance to 



M31 (~ 400 kpc) The baryonic fraction of the halo mass density should 
be about equal to the universal value, Jb — 0.04 — 0.1, or smaller, if a larger 
fraction of the baryons dissipated their energy and reside in the disk. 

(2) Suppose further that a substantial fraction of halo baryons undergo 
moderate dissipation (though not enough to collapse to form a disk), so that 
the baryonic halo remains roughly spherical but shrinks in size. Specifically, 
assume that both cdm and baryons exist in isothermal halos truncated at 
different radii, Rb and -Rcdm = PRb- If niost of the baryons are in the halo, 
the ratio of the halo mass in baryons to that in cdm is — /b) — Ib, 
and interior to Rb the ratio of baryons to cdm is PfB/i^ — fs) — Pfs- That 
is, the local baryon to cdm ratio is increased by the ratio of their truncation 
radii. To be concrete, if -Rcdm ~ 200 kpc and Rb ~ 30 kpc, then interior to 
30 kpc the ratio of baryons to cdm is about 7/s ~ 0.3 — 0.7. 

(3) Consider a very radical scenario, one where our galaxy formed by non- 
gravitational forces (e.g., hydrodynamical shock waves), so that the baryonic 
part of the galaxy was assembled first and the cdm halo accreted subse- 



quently. By using the spherical accretion model |T2| one can estimate the 
mass of cdm halo that is added. Suppose at time tt a point mass Mq is 
placed at the origin of an otherwise smooth critical Universe comprised of 
cold dark matter. Thereafter, all cdm particles in the Universe are gravita- 
tionally bound to Mq and ultimately cease moving away (clearly the model 
is only applicable within the sphere of influence of Mq, say half way to M31.) 

According to this model, cdm is accreted in a self-similar way, with a 
density profile p oc r~^/^ for r <C r^^; at time t the "turn-around" radius 
r,(t) = 0.771(GMo)i/3t8/7t2/9. The mass accreted by time t and interior to 

radius r < r, is: M(r,t) = 1.39{MorY/^/G^/Hl^\ Taking Mq ~ lO^^M© - 
(baryonic mass of our galaxy), the turn-around radius today is around 1 Mpc, 
and the mass accreted within 100 kpc is about 7 x 1O^'^M0(1 + Zi^^'^h^^'^ ~ 
IO^^Mq for ~ 1 — 2 and h ~ 0.5 {zi is the red shift corresponding to 
time ti). The density of cdm particles 8.5 kpc from the galactic center is 
0.7 X 10~25gcm-3(i + z.)3/4/ii/2 _ 10-25 gcm-3 for 2i ~ 1 - 2 and h ~ 0.5. 
Even in this radical model, the amount of cold dark matter accreted within 
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100 kpc is about equal to the baryonic mass, and the local density of cold 
dark matter is only slightly lower than the estimates in scenarios (1) and (2). 

(One can consider a more extreme version of this scenario: suppose that 
the Milky Way resides in an underdense region of the Universe, which, in 
the absence of the point mass Mq, behaves like a small portion of an Qq < 1 
Universe. In this case, the total cdm mass accreted is (1 + Zj)f2oMo/(l — f^o); 
again comparable to Mq, and so the previous comments apply.) 

Microlensing. Now we turn to what one can hope to learn experimentally 
from the study of microlensing of stars in the LMC by dark stars in the 
halo of our galaxy. Such experiments can only "detect" halo baryons if they 
exist in the form of 10~^Mq to lOOM© dark stars. Little is known about 
the kind of objects halo baryons would form (in part, because of our poor 
understanding of star formation in general). The strongest statements that 
can be made concern in what form halo baryons cannot exist and lead to the 
suggestion that halo baryons are likely to be stars of mass 10~^ Mq to O.lMg 







[0, a range that can be probed by microlensing. 

In modeling a two-component halo we assume that: (i) baryons and cdm 
exist in separate, spherically-symmetric isothermal halos, with core radii as 
and acdm and density profiles, 



P^{r) = PlocaM -2——J ' (1) 




where piocai,j is the local density; (ii) the core radii are between 2 kpc and 
8 kpc; (iii) halo baryons are dark stars of mass M. Both the baryonic and 
cdm halos play a role in determining the galactic rotation curve, but of course 
only the baryonic halo determines the microlensing rate. 

Next we compute the microlensing rate for stars in the Large Magellanic 
Cloud (LMC) as a function of the assumed local density of cdm for a suite of 
reasonable halo models (see below). That rate depends upon the distribution 
of baryonic matter in the halo and is given by 



rxmax dxyx{l — x) 

r(Plocal,cdm) Ctj) = UJqUt / . — — ; ^, (2) 



A' + Bx + x 



where Wo = V |]i^^J^^'Piocai,B, A' = {al+rD/L"^, andB = -2(ro/L) cos6cos/, 
vh is the halo-velocity dispersion, and L ~ 50 kpc, b = —33°, and / = 281° 
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are respectively the distance, and galactic latitude and longitude of the LMC 
The quantity ut, the threshold impact parameter in units of the Einstein ra- 
dius, is set by the minimum amplification that can be detected; e.g., ut = I 
corresponds to a amplification threshold of 1.34 which is typical of current 
searches 0] . The quantity Xmax is the lesser of 1 and the distance to the edge 
of the baryonic halo along the line of sight to the LMC in units of L. 

We compare all microlensing rates to a fiducial model, a baryons-only, 
"best fit" halo model with as = 5kpc, vh = 270kms~^, normalized to have 
rotation velocity of 220kms~^ at our position. The microlensing rate for this 
model is Tq = 1.66 x 10~^ut / \J M / Mq events yr^^. For further discussion of 



microlensing we refer the reader to Ref. P |. 

Our suite of models was constructed as follows: For each assumed value 
of the local cdm density, we allow the core radii to vary separately between 
2 kpc and 8 kpc; the value of piocai.B is determined by constraining the rotation 
velocity at our position to be vdro) = 220 ± lOkms"^, 



Plocal,B 



(aB/ro)tan ^(ro/as) 



X 



V 



2 



AuG 



Plocal,cdm(acdm + ^o)[l ^ («cdm/'^o) taU (ro/flcdm)] 



where v{rQ) ~ 130 ± 17kms"^ is the portion of the local rotational velocity 
that is supported by the halo. To ensure that a given model is "reasonable" 
we construct the rotation curve; in so doing we also take into account the 
contributions of the disk, bulge, and central components of the galaxy by 
using results from Ref. 0. For our "criterion of reasonableness" we follow 



Ref. |]15[: the relative difference of the maximum and minimum rotation 
velocities over the interval 3 kpc to 18 kpc must be less than 14% (most of 
the models pass this test; see Fig. 2). In the limit of a single halo component, 
the local halo density in our models is 3 — 7 x lO^^^gcm"^, consistent with 
previous estimates 0. 

The range of microlensing rates for our suite of reasonable models is 
shown in Fig. 1 as a function of the local cdm density. The microlensing rate 
is relatively insensitive to the local cdm density — the range of cdm densities 
consistent with a given F spans about a factor of 3 — and is relatively sensitive 
to the halo model parameters — for fixed cdm density F varies by around 
±50%. This is not surprising; first, most of the microlensing is due to objects 
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between 10 kpc and 30 kpc from the galactic center, so only this part of the 
baryonic halo is probed. Second, the galactic halo is not well constrained by 
rotation-curve data, so the "phase space" of reasonable models is large |TB| . 

Clearly, microlensing can only provide limited information about the local 
cdm density. For example, the microlensing rate can be as large as its value in 
the best fit, baryons-only (fiducial) model in a two-component model where 
the local cdm density is 3 x 10~^^gcm~^. Or, suppose that the microlensing 
rate were determined to be half the fiducial value; for our suite of halo models 
the local cdm density is 1 — 5 x 10~^^gcm~^. On the other hand, if the 
observed microlensing rate is found to be small, say 10% or less of the fiducial 
value, based on our set of models one could argue that the local cdm mass 
density is at least 2 x 10~^^gcm~^. 

To be specific about the dependence of the microlensing rate on the halo 
model, for a given cdm density the higher rates occur in models with larger 
baryonic core radii and larger values of the local rotation speed. In models 
with a truncated baryonic halo (not shown) F is insensitive to the truncation 
radius provided that it is greater than about 30 kpc; this is because most 
of the microlensing is due to halo objects between 10 kpc and 30 kpc from 
the galactic center. There is a small dependence upon the distance to the 
LMC which is not shown; changing the LMC distance by ±10% changes F 
by about ±5%. We should emphasize that our two-component models are 
very simple; one could easily imagine more complicated models, e.g., where 
the halos are not spherically symmetric. This increases further the range of 



plausible microlensing rate for a given cdm density |T7 



Finally, a fine point, in computing the microlensing rate, we have followed 



Ref. IIJ] in assuming that the distribution of halo velocities is Maxwellian, 
which leads to the factor of vh in loq, cf. Eq. (||). This is only strictly 
true for a galaxy model consisting solely of an untruncated, zero-core radius 
halo. In that regard, our models (like most) are not self-consistent. To 
explore the sensitivity of our results to this inconsistency, we replaced vh by 



3/2^0110(20 kpc); that is, we used the circular velocity at the radius where 
most of the microlensing occurs to characterize the halo velocity dispersion. 
This increased the sensitivity of the microlensing rate to halo model, though 
very slightly (by a few per cent). Further, we also computed the optical 
depth for microlensing, cf. Eq. (4) in Ref. which does not depend upon 



the distribution of halo velocities, and it exhibits a similar variation for fixed 
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cdm density as T. Both of these facts suggest that our estimates of the model 
uncertainties in the microlensing rates — which is our main concern — have not 
been affected significantly by lack of self consistency. Absolute rates will of 



course depend more strongly on consistency 17 



Summary. If the density of the Universe is significantly greater than 
about 10% of critical density, as a number of observations indicate and several 
arguments suggest f^, Q , then there are two dark matter problems, baryonic 
and nonbaryonic. If the nonbaryonic dark matter is cold dark matter, which 
seems to be the most promising possibility, then the halos of spiral galaxies 
should contain both baryons and cdm particles. Provided the formation 
of the halo involved only gravity, the local baryonic fraction of the halo 
should be small, less than about 10%. If baryons in the halo underwent 
some dissipation, so that the baryonic halo contracted relative to the cdm 
halo, the local baryonic fraction is increased by the contraction factor and 
could be significantly higher. In the extreme, if the baryonic portion of 
the galaxy formed first through nongravitational forces and a cdm halo is 
accreted thereafter, the mass of cdm in the halo is comparable to that in 
baryons and the local cdm density is still 10~^^ g cm~^ or so. 

Because of uncertainties inherent in modeling our halo and the fact that 
microlensing only probes the part of the halo between 10 kpc and 30 kpc from 
the center of the galaxy, it is difficult at present to learn much about the cdm 
content of own halo from microlensing. For example, if the microlensing rate 
were found to be equal to that expected in the best fit, baryons-only halo 
model, the local cdm density could be as large as 3 x 10~^^ g cm~^ when the 
uncertainties in halo models are taken into account. On the other hand, if 
the microlensing rate were found to be small, say 20% or less of the baryons- 
only model, based on our models one could argue that this is evidence for 
a local cold dark matter density larger than about 2 x 10~^^gcm~'^. While 
the MACHO and EROS collaborations have yet to discuss the microlensing 
rate that can be inferred from their candidate events [0, a naive analysis of 
their data suggests that the rate could be as low as 20% of the baryons only 
model HTH 
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1 Figure Captions 



Figure 1: The microlensing rate T for our suite of reasonable, two-component 
halo models as a function of the local cdm density, normalized to a best fit, 
one-component (baryons) halo model with core radius of 5kpc. The heavy 
curves show the extreme range of the models; from top to bottom the lines 
correspond to models with qb — Ocdm = 6.5, 5, 3kpc. 

Figure 2: A sample of galactic rotation curves for our models. The heavy 

curves correspond to extreme models (top: ttcdm = 2kpc, = 8kpc, and 
Piocai.cdm = 0.8 X lO^^^gcm^^; bottom: flcdm = ciB = 2kpc). The two middle 
curves correspond to acdm = clb = 6.5 kpc (upper), 5kpc (lower). 
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